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We have studied the iso-height lines on the WO3 surface as a physical candidate for conformally 
invariant curves. We have shown that these lines are conformally invariant with the same statistics 
of domain walls in the critical Ising model. They belong to the family of conformal invariant curves 
called Schramm- Loewner evolution (or SLE K ), with diffusivity of /t ~ 3. This can be regarded as 
the first experimental observation of SLE curves. We have also argued that Ballistic Deposition 
(BD) can serve as a growth model giving rise to contours with similar statistics at large scales. 

PACS numbers: 68.37.-d, 05.45.Df 



The study of random surfaces, especially their statis- 
tical properties as well as growth and evolution dynam- 
ics, has been growing over the last two decades. They 
describe many important problems of real surfaces ap- 
pearing in condensed matter physics such as deposited 
metal. In addition, this kind of problems is closely re- 
lated to other problems in physics such as fractured sur- 
faces, string theory, phase transitions in two dimensions 
etc 0, Q . The aim of this paper is to show that confor- 
mally invariant curves appear on the iso-height contours 
of deposited films of WO3 . This is the first observation 
of SLE in real physical system. We believe that the same 
result may be observed for other surfaces. 

Metal oxides are a large family of materials possessing 
various interesting properties. One of the most interest- 
ing metal oxides is tungsten oxide, WO3, which has been 
investigated extensively because of its distinctive applica- 
tions such as electrochromic 0, 0, H, @, H, p hotochromic 
[|| , gas sensor 0, E| , photo-catalvstfia] , and photo- 
luminescence properties [13j. 

Many properties of tungsten oxide are related to its sur- 
face structure (e.g. porosity, surface-to-volume ratio) 
and surface morphology and statistics such as grain size 
and height distribution of the sample. These properties 
are also affected by conditions during the growth process 
such as deposition method. One can change the statis- 
tics of the growth surface by imposing external parame- 
ters such as annealing temperature which can even cause 
phase transition in the sample 

On the other hand, calculation of various geometrical 
exponents of a random Gaussian surface is of interest to 
theoretical physics. From this point of view, we want to 
study the morphology and geometrical statistics of ex- 
perimental grown surfaces of WO3 . We consider contour 
lines (the nonintersecting iso-hight lines) on the WO3 
samples and show that they are conformally invariant. 
This analysis when applied to the contour lines of the bal- 
listic disposition (BD) model shows similar conformally 
invariant surface in large scales. The scaling behavior 
of contour lines in the growth models are widely investi- 
gated. For example in [15J the authors derived some of 




the universal exponents of contour lines, especially their 



FIG. 1: (color online). AFM image of "as deposited" sample 
of WO3 thin film in a scale of lfim x 1 /j,m with the resolution 
of 1/256 fjm. 



fractal dimension Dq and its relation to the roughness 
exponent a. Moreover, Schramm and Sheffield recently 
showed that the contour lines in a two-dimensional dis- 
crete Gaussian free field are conformally invariant and be- 
long to a wide class of conformally invariant curves called 
Schramm-Loewner evolution (SLE K ) curves, whose dif- 
fusivity k converges to 4 [l6j- every conformal curve in 
two dimensional plane can be studied by SLE [l7| (to re- 
view about SLE see [111). Moreover, it has been proved 
by Smirnov that some domain walls in the critical Ising 
model, in the scaling limit, can be described by SLE [l9| . 
In addition some physical systems have been studied re- 
cently using some numerical methods. It has been shown 
that the statistics of zero-vorticity lines in the inverse cas- 
cade of 2D Navier-Stokes turbulence displays conformal 
invariance and can be described by SLEq which is in the 
universality class of percolation |20| . Similar studies have 
been done for inverse cascade of surface quasigeostrophic 
turbulence [2l[ , domain walls of spin glasses [i^] and the 
nodal lines of random wave functions [23[ . 

To study the geometrical statistics of the surface of 
WO3 , 34 samples were independently deposited on glass 
microscope slides with the area 2.5cm x 2.5cm using a 
thermal evaporation method in same conditions. The 





FIG. 2: (color online). The connected domains with posi- 
tive heights on the WO3 surface. Negative height regions are 
black. 



FIG. 3: (color online). Some of the positive height cluster 
boundaries on a WO3 surface (see also Fig. 0. 



deposition system was evacuated to a base pressure of 
~ 4 x 1CP 3 Pa. The thickness of the deposited films 
was considered to be about 200 nm. measured by the 
stylus and optical techniques. Using the atomic force 
microscopy (AFM) technics we have obtained the height 
profiles h(r) of the rough surfaces with the resolution of 
1/256 \im in the scale of 1/im x 1/im. The AFM scans in 
this scale were performed in the various non-overlapping 
domains (10 profiles for each sample) from the centric 
region of the deposited samples (Fig. [TJ-to ensure that 
the boundary effects are negligible for the AFM profiles. 
First we have determined the clusters by considering the 
connected domains with like-sign heights (mean height is 
set to zero for a range of heights around the mean value 
much less than h rms )(Fig. [2]), then the cluster boundaries 
(contour lines) have been specified (Fig. [3]). We have ob- 
tained an ensemble of such contour lines using MATLAB 
algorithm for contour plot. The calculation of the fractal 
dimension of these contour lines yields Do = 1.38 ± 0.02 
using the box-counting method. The roughness exponent 
for the samples is calculated by using the second-order 
structure function analysis, which gives a = 0.33 ± 0.03. 

In order to investigate the scale dependency of the 
contours we have analyzed their multifractal behavior. 
The generalized fractal dimensions are defined by D q = 
pjfim^o!^^ where Z q (a) = Pt{a) q , and the sum 
runs over all the non-overlapping boxes with size a which 
cover the curve. Pi(a) is the mass of the curve in the zth 
box. The multifractal spectrum of the contours is inde- 
pendent of q with the value of 1.38 ± 0.03. It shows that 
the curves are scale invariant. To investigate more the 
scale dependency of contour lines, one can check some of 
the exponents associated with clusters and loops. Fig. [4^ 
shows the scaling of the average area of clusters M with 
their radius of gyration R, M ~ R D which yields the frac- 
tal dimmension of clusters D = 2 — (8 — k)(3k — 8)/32k 
24 1 . Fractal dimension of loops relates their length s to 
radius of gyration R as s ~ R D ° . In our case, it is in good 



agreement with both the fractal dimension obtained us- 
ing the box-counting method and that of domain walls 
of critical Ising model (Fig. [4Jd). We have also checked 
the other scaling relations, the average number of loops 
with radius R, n{R) ~ R- 2+a (Fig. Bfc) and that with 
perimeter s, n(s) ~ s (- 2 + a )/ D o where a is the rough- 
ness exponent (Fig. QJi). Another is the scaling relation 
for the number of loops with the area greater than A, 
n(A) ~ A- 1+a ' 2 . 

Consistency of the results for the contour lines of WO3 
samples and Ising cluster boundaries suggests that a con- 
tour line of these two belongs to the same universal- 
ity class. However assumption of translational and ro- 
tational symmetry within contour ensembles has been 
made, which seems right. Therefore the iso-hcight con- 
tours of WO3 may be conformally invariant and belong 
to the SLE K curves with k = 3. The fractal dimension 
Da = 1 + k/8 is also consistent with our finding provided 
k = 3. 

To examine this hypothesis directly, we first consider 
an arbitrary placed horizontal line representing the real 
axis in the complex plane across the AFM profiles. Then 
we cut the portion of each curve 74 above the real line as 
it is in the upper half plane H. 

Now let us to consider the chordal SLE in the H 
from to 00 which is described by the Loewner equa- 
tion dtgt(z) = gt (,)_£ t ■ If we parameterize the contour 
lines in the H with the dimcnsionless parameter t, gt{z) 
maps the upper half plane except the curve up to time 
t, to H itself. For any self avoiding curve 74 the driving 
function £t is a continuous real function with the corre- 
spondence of g t (j t {z)) = £t- In the other words, gt{z) 
in each time, maps the tip of the curve to a real point 
£ t . Random curves whose driving functions are propor- 
tional to the Brownian walk B t are conformally invariant. 
So, to investigate the conformal invariancy of the contour 
lines we should explicitly show that their driving function 
Ct = s/ftBt i.e., its increments are identically and inde- 
pendently distributed and (£ f 2 ) = nt, where k is diffusion 
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FIG. 4: (color online). Cluster and loop statistics for iso- 
height lines of WO3 samples in scales 1/im x 1 /im (Squares) 
and 5/im x5 /im, lO^tm xlO /mi and 20/mi x20 /mi. The size 
of the samples are rescaled to be 1 with 256 2 lattice points, (a) 
The average area M of clusters versus the radius of gyration 
R. (b) The perimeter of loops s versus the radius of gyration 
R. (c) The average number of loops with the radius of R. 
(d) The average number of loops with the perimeter s. (e) 
Number of loops of area greater than A. The error bars are 
almost the same size as the symbols in the scaling regions, 
and have not been drawn. These results are very close to the 
corresponding statistics of domain walls in the critical Ising 
model (the solid lines). Note that in Figs, (c), (d) and (e) the 
change in the exponents relative to the critical Ising model is 
due to the roughness exponent [l5( |. 



constant. 

There are many numerical algorithms for extracting 
the driving function of a given random curve by inver- 
sion of loewner equation [25j . We have used the succes- 
sive discrete, conformal slit maps - based on the piece- 
wise constant approximation of the driving function - 
that swallow one segment of the curve at each time step 
p3l ]. To do this at first, we have determined all of the 
contour lines 7tS of the WO3 samples as described above 
by sequences of the points {zq, Z\, ■ ■ •, zjv}- we obtained 
2319 such curves with the average number of points about 
1210. After that all of the contour lines were mapped by 
ip(z) = znz/(zn — z). To avoid numerical errors only 
the part of the curves corresponding to capacity t < 0.25 
were used. 

Then, by setting the starting point zq in (0,0), all of 
the points except the first one have been mapped with 
g t {z) = V^-CO 2 + 4i + 6 with t = \lm{ Zl ) 2 and 
£t = Re(zi) and the resulting points have been renum- 



FIG. 5: (color online). Statistics of the driving function £(t). 
main frame: the linear behaviour of (£(t) 2 ) with the slope of 
k = 3±0.2. Lower-right inset: The probability density of £(t) 
rescaled by its variance k at different times t = 0.04, 0.06, 0.08 
is Gaussian. Upper- left inset: The correlation function of the 
increments of the driving function 8^(t). 

bered by one element shorter. The algorithm is repeated 
until the whole curve is transformed. The final output of 
the process is a series of driving functions obtained from 
the contour ensemble of the WO3 samples. (The accu- 
racy of the code has been checked on an ensemble of SLE 
traces with known k where it yielded the correct value 
with an error of ~ 3%). 

As shown in Fig. the statistics of the driving function 
-considering finite size effects- converges to a Gaussian 
process with variance (£(£) 2 ) = K t an d « = 3±0.2 and mi- 
correlated increments. It confirms that the contour lines 
of the WO3 samples are locally SLE3 curves i.e., con- 
formally invariant. To investigate the effect of the scale 
of AFM images on the statistics of the contour lines, we 
took over 10 2 AFM profiles from samples at each of three 
scales 5/imx5 /im, 10/imxl0 /mi and 20/mix20 /jm with 
a fixed number of lattice points 256 2 . After rescaling of 
size to one, we observe no conclusive differences within 
numerical errors (Fig. [4j. 

In the following, due to some similarities between some 
of the results for WO3 surface and BD model, we wish 
to briefly discuss this model. 

To investigate the behavior of the contour lines of BD 
growth surface, we simulated 10 2 independent samples 
with the size of 1024 x 1024. The roughness exponent 
of the samples is 0.36 ± 0.02 [II]. As shown in Fig. [6ja,, 
the statistics of the heights are self-similar and the height 
increments have approximately Gaussian statistics [261 ]. 
The multifractal spectrum of the BD contour lines has 
been numerically obtained which shows that the contour 
lines seem to be scale invariant. The fractal dimension of 
the lines is 1.39 ± 0.03 (Fig. [6)3). In the first view, these 
results show that the BD's contour lines may belong sta- 
tistically to the same universality class as the domain 
walls in the critical Ising model. Applying the Inverse 
loewner equation, we have analyzed the driving function 
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tics of the contour lines are conformally invariant and it 
is very similar to the statistics of the domain walls in the 
critical Ising model. However, this approach can be used 
for other experimental samples to know more about the 
relation of the statistics of the contour lines to the var- 
ious geometrical and physical properties of the surfaces. 
Moreover, we have proposed a model of growth , ballistic 
deposition model whose contour lines seem to belong to 
the same universality class as the Ising model at large 
scales. 

We wish to thank A.Z. Moshfegh, R. Azimirad and S. 
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thank B. Fallah Hassanabadi, S.M. Fazcli, R. Scpchrinia 
and especially R. Nourafkan and M. D. Niry for their 
useful discussions. We also wish to thank F. D. A. Aarao 
Reis, J. Cardy, O. Schramm and R. Stinchcombe for pro- 
viding helpful comments on this manuscript. 



FIG. 6: (color online), (a) Second-order structure function of 
height differences and, in the inset, rescaled and normalized 
probability density functions for r=0.02, 0.04 and 0.06 com- 
pared to a Gaussian density (solid line). Data have been ob- 
tained from a simulated BD growth surface in square domain 
of size 1 with 1024 2 lattice points, (b) An instant contour line 
and, in the upper-inset, generalized multifractal dimensions 
for the BD contour lines which are independent of q within 
statistical errors. The dashed line is set in the 11/8. Lower- 
inset shows the scaling of the perimeter of the loops s versus 
the radius of gyration R compared with the one in the Ising 
model (dashed line). 



of the iso-height curves. This analysis shows that in some 
large scales the increments of the driving function are in- 
dependently distributed with Gaussian statistics of vari- 
ance k = 3.1 ± 0.2. 

A point which remains is the difference between the 
standard way of generating an interface in the Ising 
model and the one considered in this paper. To gen- 
erate chordal SLE in the half plane for the Ising model, 
one usually takes fixed boundary conditions on either side 
of the origin. Whereas in the procedure proposed here, 
we choose an arbitrary horizontal line which clearly has 
more complex boundary conditions on it. Certainly dif- 
ferent boundary conditions can introduce different drift 
terms in the driving function. Question is which bound- 
ary condition with the above procedure leads to chordal 
SLE in the Ising model. We have run simulations on 
the two dimensional critical Ising model with periodic 
boundary conditions on the square lattice of size 1024 2 , 
and found that the resulting driving function has leading 
term (considering finite size effects) with the same statis- 
tics of Brownian motion with diffusivity of k = 2.9±0.2. 
Details will appear in a forthcoming work. 

In conclusion, we have found an experimental exam- 
ple, the iso-height lines of the WO3 surface, for SLE 
curves. It has been numerically shown that the statis- 
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